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Abstract 

We analyze the dynamics of a network of electrically coupled, modified 
FitzHugh-Nagumo (FHN) oscillators. The network building-block architec¬ 
ture is a bidimensional squared array shaped as a torus, with unidirectional 
nearest neighbor coupling in both directions. Linear approximation about 
the origin of a single torus, reveals that the array is able to oscillate via a 
Hopf bifurcation, controlled by the interneuronal coupling constants. Group 
theoretic analysis of the dynamics of one torus leads to discrete rotating 
waves moving diagonally in the squared array under the influence of the 
direct product group x x 7^2 x E>^. Then, we studied the existence 
multifrequency patterns of oscillations, in networks formed by two coupled 
tori. We showed that when acting on the traveling waves, this group leaves 
them unchanged, while when it acts on the in-phase oscillations, they are 
shifted in time by 0. We therefore proved the possibility of a pattern of oscil¬ 
lations in which one torus produces traveling waves of constant phase shift, 
while the second torus shows synchronous in-phase oscillations, at N— times 
the frequency shown by the traveling waves. 

Keywords: group-theoretic, coupled-tori, traveling-waves, in-phase 
oscillations, Hopf bifurcation. 
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1. Introduction 

The phenomenon of multifrequency oscillatory patterns has been initially 
described in electrical systems, where one or many variables happened to 
oscillate with two or three times the frequency of the other oscillators when 
there was no obvious symmetry reason to do so Q (and references therein). 
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This phenomenon received a solid groim theoretic treatment by the works 

3 ,( 3 , 0 , 0 , 0,0 


of Golubitsky and Stewart 

Chossat in [l2| and it has also been conhrmed experimentally in 
and 


Concretely, Palacios et ah in [17[ , [1^ proved the existence of such 


17 , 16 


multifrequency patterns in electrical systems of Duffing oscillators. 

Duffing oscillators where the first dynamical systems used to model the 
brain activity and in particular, to reproduce the electroencephalograms (3 - 
However, the first mathematical model for the electrical signaling or firing for 
individual as well as coupled neurons, developed in direct relationship with 
the underlying neuronal physiology, was the Hudgkin-Huxley model [26| and 
its simplified version, the FHN model (3, (3. Neurophysiological informa¬ 
tions on brain processing indicate that neural networks must function in a 
relatively narrow frequency range 0. This contrasts with the observation 
that brain uses a wide frequency spectrum including several oscillation ranges 
(delta, theta, alpha, beta, gamma), which gives poor efficiencies of signal cod¬ 
ing by different frequencies, because of low signal resolution in this frequency 
range 0. It has been proposed that hierarchical feature integration can be 
accomplished using multifrequency quasiperiodic oscillations 0. For a two- 
level system, the authors propose that synchronization of high-frequency os¬ 
cillations is used to integrate simple features 0, (3) whilst for more complex 
features out-of-phase oscillations of constant phase shift were proposed ^ 
and 3 . Multifrequency oscillatory patterns are in fact widely described i 
neural systems. Many periodic activities are ori gin ated b 
forming central pattern generators 


m 


roups of neurons 
11| . In addition, co¬ 
existence between in-phase oscillations and trave ling waves of constant phase 


20 


a, 10 


ling 

shift have been reported in the nervous system 3 references therein. 

Symmetry is a constantly present feature in nervous system architecture. 
The already mentioned central pattern generators, visual cortex neuronal 
circuits in mammals 


30| . [23[ or plasmodial slime mold [2^ are just a few 


examples in which symmetry groups act on the dynamical systems repre¬ 
sented by biological oscillators 0 . 0 and 0 . Consider for instance the 

following system of differential equations, — = f(x) where f : M”' —)■ M"’ is 

at 

a smooth vector held. Let F be a group that acts in M"". Let us denote the 
action of F on the space of vector helds on M” by F also. In general we have 
the result that if the action of a symmetry group F on an ordinary differential 

G?X 

equation — = f(x,/i), equivariant under the action of the symmetry group. 
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where x G M"", /r G M and f is smooth, then 7f(x, fi) = f (7X, /i), V 7 G F. From 
the nniqneness of the solutions of initial value problems, it follows that either 
x{t) and 7 x{t) are disjoint trajectories resulting from a new periodic solution 
or they differ from each other by a phase shift. This leads to the following 
conclusion: when a system inherits a symmetry and when it undergoes Hopf 
bifurcatiom resulting periodic solutions also inherit certain symmetry prop¬ 
erties [H, jH- We study in this paper, the dynamics of a symmetric network 
of coupled neurons, that is, a network on which act certain symmetry groups. 
We focus our analysis on the multifrequency patterns of oscillations which 
are dependent on the network architecture. 

In our approach, the neurons are of FHN type. For a single neuron, the 
FHN model [2^ is modified according to [l^; the variables assumed in this 
model are the membrane potential x and (a surrogate for) the ionic current 
y. The state of the neuron is thus specified by a point {x,y) G and the 
internal dynamics is 

X = ax — x^ — y 

■ h 1 

y = ox — cy 


where a, b are parameters, and 0<a<l,6>0,c>0 and 6 > a. In a two-cell 
network, the dynamics of one cell is modified by coupling influences from the 
other cell. This can be imagined, as in 2l|, by adding an applied current 
to the X equation, this being a function of the state of the other cell. Thus 
xi = {vi,wi),X 2 = {v 2 ,W 2 ), and g{u,v) = {axi-xl-yi + k{xi-x 2 ),bxi-cyi), 
where k represents coupling strength. 

The paper is organized as follows. In section |2] we describe the building 
block structure of our network: a three-dimensional lattice shaped as a torus, 
and identify the symmetry group acting on the coupled differential systems 
located at the nodes of the lattice. In section |3l we start the analysis of 
the network formed by one torus. We analyze the dynamics of the linear 
approximation of system (| 2 ]) about the origin, by means of the explicit ex¬ 
pressions for eigenvalues and eigenvectors, and identify the existence of limit 
cycles arose from the Hopf bifurcation which depends on the interneuronal 
couplings. In section 0] we continue to study the dynamics of one torus, this 
time from a group theoretic angle. We apply the Theorem [1] from [1^, to 
identify the two-dimensional isotropy subgroups of the symmetry group act¬ 
ing on the global system |2l and prove the existence of N— distinct discrete 
rotating waves which are diagonally translated in the lattice. In section 01 
we carry our analysis of the dynamics of two coupled tori. We show that 
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certain groups of symmetries can be associated with periodic patterns where 
an entire array oscillates at different frequencies, and apply this argument to 
describe the multifrequency patterns of oscillations in two coupled tori. 


2. Network structure 


Palacios et ah as well as Longhini et al. 0 experimentally proved the 
existence of multifrequencies in unidimensional rings of overdamped Duffing 
oscillators as well as Van der Pol oscillators 0- Since our main aim here 
is studying the patterns of multifrequencies, it appears natural not only to 
give a detailed description of our network structurejbut also to briefly com¬ 


pare it with the cited network architectures. In and [17[, it has been 


analyzed the dynamics of two unidimensional and oriented rings of cells with 


nearest-neighbor intrarray coupling and all-to-all interarray coupling. In [16 


they considered arrays formed by N oriented rings of cells, the network ar¬ 
chitecture respecting unidirectional coupling in each array and also between 
adjacent cells in the N arrays. 

In particular, the general network presented in [l^, is such that the hrst 
row (or column) is not connected to the last row (or column). Their net¬ 
work is constructed from a rectangular array after the top and the bottom 
edges are identihed with each other, such that the hnal structure respects 
the cylindrical topology. They used opposite directions of the intraarray 
couplings. 

In this work, we start to build our network, by considering a squared 
lattice, the nodes of which represent FHN neurons, N in each direction. We 
hrst fold the squared array into the shape of a cylinder, followed by closing 
up the ends of the cylinder. Thus, the building block of our network is a 
two-dimensional array wrapped in the shape of an empty torus. A similar 
structure served Miller et al. [l^ to model an array of Van der Pol oscillators 
as the cartesian product of three rings. Our interarray couplings have same 
directions; we do not use cross couplings. 

As a hrst approximation we assume that all neurons are identical. Each 
neuron of the network is coupled to the two forehead neurons (of its four 
nearest neighbors), one in each direction, and we respect the same direction 
of interring couplings. For example, neuron located at («,/S) in the array, is 
coupled to neuron {a + 1, (3) along the a direction, and with neuron (a, (3 + 1) 
along the 13 direction (see Figure H]). The dynamics of the neuron at the 
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location is described by the variables Xc^i 3 ^ya ,/3 


Xa,l3 ^a.,13 “t" T^^(^a+l,/3: “1“ 

ya.,13 ^ya,0 


( 2 ) 


where 

^(^a+1,/9) ^a+1,0 “t“ ^a,0i /q\ 

h(^a,/3+l) 2^0!,(9+1 Xa,Pi 

are the interneuronal intratorus coupling functions ,7 is the coupling constant 
in a direction, while <5 is the coupling constant in 13 direction. 



unidirectional intra-ring 
nearest-neighbor coupling 
on torus surface 



(b) 


Figure 1: (a) Schematic diagram of the bidimensional array on the torus surface: along 
the a direction (the coupling constant is 7 ), neuron (a + l,/3) is unidirectionally coupled 
with neuron {a,/3), neuron (a,/3) is coupled with neuron {a — l,/3), etc. Along the (3 
direction (the coupling constant is (5), neuron (a,/3 + 1) is unidirectionally coupled with 
neuron (a, /3), neuron (a, j3) is coupled with neuron (a, /3 — 1), etc.; (b) Schematic diagram 
of a 2 -tori network obtained by bidirectional inter-tori coupling (the coupling constant is 
e). 


In the following, we will try to identify and describe the symmetry group 
of one torus. Let us denote by Zjy, the group of cyclic permutations of N 
neurons in each direction on the torus surface. In this torus-like neuronal 
network with oriented nearest neighbor coupling, the symmetry group is not 
Zat X Zat but Zat X Zat X Z2, where the additional symmetry is introduced by 
the property {x,y) 1 —)■ {—x,—y). More specihcally, the action of this group 
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on 0 with coordinates {x,y) is: 

(t{xi^I3, . . . , Xn,i3] 2/1,/3, • • • , yN,l3) = (a^2,/3, • • • , Xn,i3, 3^1,/?; 1/2,0, • • • , yN, 0 , 2/1,/?) 

• • • ) ^a,N) ya,l, ■ ■ ■ ■> ya,N) {,-^a,2, • • • , ^a,N, ^a,l, 2/a,2) ■ ■ ■ , ya,N, ya,l) 

^{Xl, 0 , • • • , XN^j 3 ] 2 / 1 ,/3, . . . , yN, 0 ] Xa,l, ■ ■ ■ , Xa,N] 2/q, 1) • • • ) 2/a,Ar) = 

( ^l,/3) • • • ) X]\f^ 0 , 2/1,/9) • • • ) yN, 0 ) Xa,l, ■ ■ ■ 1 Xa,N-I ya,!, ■ ■ ■ ■, ya,N) 

where = (a), Z^ = (p), and Z 2 = (ru). Let F = Z^ x Z^ x Z 2 . 

In a network with cylindrical topology, coupling functions k {xa+i, 0 ,Xa, 0 ) 
and /i{xa, 0 +i,Xa, 0 ) fail to be F equivariant for the reason that neurons in 
the first row (or column) are not treated as neighbors of neurons in the last 
row (or column) and vice-versa. However; if the two free ends of the cylinder 
are connected to eachother, and thereby declaring first and last rows (or 
columns) as adjacent, then the resulting i<i{xa+i, 0 ,Xa, 0 ) and /i{xa, 0 +i,Xo,, 0 ) 
turn out to be F equivariant. 

In our attempt to study the multifrequency patterns in two coupled tori, 
we would like to analyze first, the dynamics of a single torus. The main 
question we try to answer to in the next section is the following: how does the 
coupling functions/parameters affect the network dynamics? Linearization 
about the origin leads us to the following. 



Figure 2: Plot of Xi^i versus xi ^2 from the integration of system corresponding 
to a torus of 3 x 3 neurons: (a) low coupling case, 7 = 0.1 and 5 = 0.1. The 
trajectory describes an open curve from the initial conditions to the stable equilib¬ 
rium point; (b) high coupling case, 7 = 2 and 6 = 2. The trajectory approaches 
a closed curve indicating the closeness of a periodic solution. The rest of the pa¬ 
rameters and initial conditions for all variables in both cases are a = 0.01, b = 
c = 0.9, x° = (0.8462,0.2026,0.8381,0.6813,0.8318,0.7095,0.3046,0.1934,0.3028),?/° = 
(0.5252,0.6721,0.0196,0.3795,0.5028,0.4289,0.1897,0.6822,0.5417). 


3. Linear analysis 

To analyze bifurcation of the system (j2]) we need to understand the eigen¬ 
values and eigenvectors of the linear approximation about the origin. Origin 
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is a singular point of ([2]) and linearization leads to the N x N block circulant 
matrix M 


M = 


A B 

0 A 
0 0 


0 

B 

A 


... 0 0 
0 ... 0 
B ... 0 


B 0 ... 0 0 A 


where A is an iV x iV block circulant matrix and 5 is an iV x iV block diagonal 
matrix. 


'd 

E 

0 


0 

0 


'f 

0 

0 


0 

0 

0 

D 

E 

0 


0 


0 

F 

0 

0 


0 

0 

0 

D 

E 


0 

B = 

0 

0 

F 

0 


0 

E 

0 


0 

0 

D 


0 

0 


0 

0 

F 


Matrices B, E and F are given by 


D = 


d 

b 



-7 

0 



-(5 0 
0 0 


where d = a + j + 6. Let us dehne G(r,s) = D + + (sF, where Cr = e n"' 

and Cr = are the and roots of unity, respectively, with 0 < r, s < 
— 1. To hnd the eigenvectors of M, we hrst start with the eigenvectors of 


where 


G(^r,s), which are 2 = 1 , a + 7 (1 - Cr) + (1 - C) - 
A(r,s)i 2 eigenvalues corresponding to each pair (r, s). Then dehne 


'^(.r,s)i 2 ~ 




1 , 2 ’ 




1,2 


To hnd the eigenvectors of M, we proceed as in [2^ so that for each pair 
(r, s), we obtain the corresponding pair of eigenvectors 


^{r,s) 


1,2 


1 

N 


^{r,s)i 2^ C ^{r,s)i 2^ C ^{r,s)i 2^ • • • ) ^ 


where 0 < k < N — 1. 
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Figure 3: Time series of variables a) Xi.i (continuous line), X 12 (dashed line) , (dotted 
and dashed line); b) X 2 ^i (continuous line), ^2,2 (dashed line) , ^2,3 (dotted and dashed 
line) and c) x^^i (continuous line), a;3_2 (dashed line) , 0:3^3 (dotted and dashed line) from 
the integration of system © corresponding to a torus of 3 x 3 neurons. Coupling constants 
are, 7 = 2 and 5 = 2, while the rest of the parameters and initial conditions are the same 
as in Figure [2] There is a constant phase shift between the variables belonging to the 
same row. Power Spectrum Density indicating that all variables oscillate with the same 
frequency. 


It is straightforward to derive the explicit expression of the 2N^ eigenval¬ 
ues of M] there is a pair of complex conjugate eigenvalues of the form 

\r,s)i 2^0 + O + 7(1 “ Cr) + < 5(1 — Cs)] 

1 /- - - ( 4 ) 

=*=2 V ® ~ ~—46. 

corresponding to each pair (r, s). 

If Ai ^2 correspond to the mode (r, s) = (0, 0) then expression (|1]) reduces 
to the eigenvalues of system ([T]) of uncoupled neurons, linearized about the 
origin, 

_ -c + a± y/(c + a)^ -46 
’ 2 

Since 7 ^ 0 and 7 ^ 0 it follows that if taking a or c as bifurcation 
parameter, the pair of complex conjugate eigenvalues ([S]) will cross the imag¬ 
inary axis of the complex plane at c = a whenever < 6 , and system ([ 2 ]) 
undergoes a Hopf bifurcation within the mode (r, s) = (0, 0). Notice that the 
Hopf bifurcation does not depend on the interneuronal coupling constants. 

We will now analyze the conditions for Hopf bifurcation of system (1^ , 
when not both r and s are zero simultaneously. Let 7ZS be the set of combi- 














nations (r, s); moreover, the set MZ = IZS \ (0, 0), defines the combinations 
set of (r, s), where not both r and s are zero simultaneously. Let us define 


■A ‘Ittv 

Oi = (a + c) + 2 (a + c) I 7 — 7 cos —— + 5 — 5 cos 


N 


2tis 

If 


and 


+ 7 

+ 5' 

+ 275 

-45, 


— 2 (fl + c) 

+ 7 

+ 275 


27rr dvrr 

1 — 2 cos —— + cos —— 

N N 

. „ 27rs dvrs 

1 — 2 cos-h cos- 

N N 

2 nr 2713 

1 — cos —— cos + cos 


27rr . 27rs 
-7 sin -5 sm 


( 6 ) 


27r (r + s) 
N 


N 


27ir Attt 
-2 sin —— + sin 


N 


sm 


N 

2'Kr 

If 


— sm ■ 


N 

2713 

Iv 


27rs dvrs 

-2 sin + sin 


N 


N 


( 7 ) 


+ sin 


27r (r + 3 ) 
N 


In the following Proposition we summarize the conditions for Hopf bifur¬ 
cation and stability of system ([2]), when (r, s) G AfZ. 

Proposition 1. Consider the linearization of system ^ around (0,..., 0), 
and assume the square root of eigenvalues in equation is a complex number 
of the form fgffh, where in general, h 7^ 0. Then, if (r, s) G AfZ, 

(i) Not both real parts ReAp^s)^^) of the eigenvalues corresponding to the 
pair (r, s), can be zero simultaneously. 




7^ 0 or 




(a) If taking as bifurcation parameter'y or 6, and y- ^ 

0 then a Hopf bifurcation occurs when —c+a+7 (l — cos +5 (l — cos 


± 


'^\l V^+ [ ) 


(Hi) Origin is a stable equilibrium of system ([^ if —c + a + 'y (l — cos 2=) + 
5 (1 — cos < min < f af + b^ + ai \ , unstable otherwise. 


27r5 \ 

N ) 
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Proof. 


(i) Let us suppose the contrary, i.e. that we have ReA(r,s)^ = ReAp,s )2 = 0. 
We will prove that this leads to contradiction. In order to evaluate the 
global real part as well as the global imaginary part of eigenvalues of 
the extended system ([2]), we need to rewrite equation (jl]) by getting rid 
of the square root. 

For this purpose, we infer a well known result from elementary algebra; 
we have that if g and h are real (h 7^ 0 ), then yjg + ih = 02 + 62 *, where 
02 and 62 are real and given by 

02 = Vg^^Th^ + g, ( 8 ) 

and 

h = (9) 

Therefore, if y/g~+Jh = \J[c + a + 7(1 — Cr) + 6(1 — C)]^ ~ 46, then 
(7 = Oi, where Oi is given in equation (j 6 ]), and 6 . = 61 is given in 
equation ([7]). 

From equations (j4]), ([6]), ([7]), (jH]) and (jO]), we have that the global real 
and imaginary parts of eigenvalues (jl]) are given by 


Re A 


(us)l,2 


[-C + O + 7 (1 - COS ^) + 6 (1 - COS ± O 2 ] 

2 


and 


ImA 


(U^*)l,2 


-7 sm ■ 


27rr 

w 


0 sm —— ± 0 - 
N 


respectively, where 02 and 62 are functions of oi and 61 given by equa¬ 
tions (IH]) and ([9]), respectively. It follows that we can rewrite equation 
dH) as 

A(r’,s)]^ 2 RsAp^s)^^ 2 + 2' 

Since we supposed that ReA(r,s)j^ = ReA(r,s )2 = 0, from equation flTOj) 
we have 02 = 0 , which means (from equation (IH])) that 61 = 0 , which 
contradicts the initial assumption that h ^ 0. 
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(ii) If applying condition ReAp,s)^ = 0 or ReA(r,s )2 = 0, to equation ([TU]) . 
we obtain the desired inequality. 

(hi) If we put the condition ReAp,s)^ ^ < 0 equation ([TU]) . we obtain —c + 
a + 7 (l - cos + (5 (l - cos < min \/a\ + h\+ . 

■ 

The linearized flow in possesses periodic solutions if ImA(r,s)^ 2 ^ 

Suppose that 7 ^ 0 or ^ 0 in ffT 0 |) . Then, Proposition [T] 

tells that by varying the coupling parameters one can qualitatively modify 
the periodic behavior of the coupled system (|2]), with respect to uncoupled 
neurons O- A remarkable property of the coupled system resides in the fact 
that self sustained oscillatory behavior of the global system can appear even 
if individual neurons taken separately act as damped oscillators, as seen in 
Figure 121 This behavior is determined by the value of the coupling constants 
7 , 5, as well as N. Let 03 = ^[ 7(1 — cos + 5 (l — cos ± 02 ] . The 
mechanism of this coupling-induced difference of the system behavior can 
be better seen if directly comparing the ReA(r^s)^ ^ from equation flTOl) of the 
coupled system ((21) with Re\i ^2 from equation ((5]) of the one-neuron-system 

o 

ReA(r,s)j^ 2 ~“ 1 “ ®3- ( 11 ) 

From equation flTTl) . we have that if ^ > 0 or ^ > 0, system loses stability 
by increasing 7 , 5, by approaching the Hopf bifurcation point. 

When N is even, however, a condition for the Hopf bifurcation of the 
discrete rotating waves, is taking into account not only nearest-neighbor but 
also next-nearest neighbor couplings 0 . 0 . 

An interesting question is determining the type (supercritical or subcriti- 
cal) of the Hopf bifurcation obtained by linearization about origin, described 
in this section. The hrst step in this analysis consists in studying the stabil¬ 
ity of the limit cycle arose from the Hopf bifurcation in the three-parameter 
family of modihed FHN systems ([I]). Hopf bifurcation of the classical FHN 
system has been studied by many authors (see for example |^), but 
these results are not applicable to our equation ([T]), which is a modihed FHN 
system. In the Appendix we apply a stability criterium due to Guckenheimer 
and Holmes [^, to prove that the Hopf bifurcation of system ([T]) obtained 
by linearization about origin, is supercritical. The stability of the limit cycle 
experimented by the coupled FHN systems within one torus (represented by 
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equation (j2])), as well as coupled tori with different intertori architectures (in 
the case of two tori, represented by equation flT5]) h is an extense analytical 
study involving reduction to the central manifold, as well as a bifurcation 
numerical analysis and it will be presented in a forthcoming paper. 


4. Analysis of the dynamic behavior of one torus with X Z^v X Z 2 
symmetry 


In a dissipative system such as ([2]), one method for finding periodic so¬ 
lutions of a given type is Hopf bifurcation. This has been analyzed in the 
previous section, from the point of view of the eigenvalues obtained at the 
linearization about origin. In the following, we would like to analyze the 
periodic solutions to system ([2]), from a group theoretic angle. A very useful 


tool is Theorem 3.4 in 14 


Theorem 1. Let F he a finite group acting on There is a periodic solu¬ 
tion to some T-eguivariant system of ODEs on with spatial symmetries 
K and spatio-temporal symmetries H if and only if 


(a) H/K is cyclic; 

(h) K is an isotropy subgroup; 

(c) dim Fix(J’f) > 2. If dim Fix(J’f) = 2, then either H = K or H = iV(iF); 

(d) H fixes a connected component ofFix{K) \ Lk- 

When these conditions hold, there exists a smooth F-equivariant vector field 
having a periodic solution with the desired symmetries. To find the kinds of 
periodic solutions in ([ 2 ]), we need therefore to determine, up to conjugacy, all 
isotropy subgroups K having dim Fix(iF) > 2 and all subgroups H for which 
H/K is cyclic. In other words, all nontrivial periodic solutions of system (12!) 
correspond to the non trivial zeros in the two-dimensional hxed point space of 
some isotropy subgroup S C ZjvxZArxZ2xS^. In the process of identifying 
these isotropy subgroups, we will first reduce F x to ITfi x Z^ x by 
considering a group homomorphism [l^; then, we will show that the only 
isotropy subgroup S with dim Fix(S) = 2, is the kernel of Z^ x x 
But before hnding these two-dimensional isotropy subgroups, let us make the 
following considerations. 
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Recall first, that we have F = x x Z 2 . If now we call Fi = Z^ x Z^, 
then we have F = Fi x Z 2 . For now, let focus our attention on the direct 
product group Fi. Let be a generator of Z^v- Its elements are written 
multiplicatively, in the form and let u) = From Z^ = (a), we have 

= 1 and from Z^ = (p), we have = 1. By following the ideas presented 
in [ 1 ^ , the two-dimensional irreducible representations am and pm of Z^ and 
Z^ respectively, over C are ■Ci) • • • > 'C[v/ 2]5 dehned as 

^z = uj"^z, m = l,2,...,^. (12) 


.^0 is the trivial representation on M. When N is even, ^isi /2 is the represen¬ 
tation on R in which both a and p act trivially. In all other cases, (where 
0 < k < N — 1), is the representation on R^ = C in which both a and p act 
as multiplication by 

The space R^^ = (x) decomposes into F 1 irreducibles according to R^^ = 
xo©,..., ©X[Ar 2 / 2 ], where the action of Fi on Xk is isomorphic to x and 
the action of zu is by —1. Similarly, the space R^^ = (y) decomposes into Fi 
irreducibles according to R-^^ = po©) • • •) ®y[N^/ 2 ], where the action of Fi on 
yk is isomorphic to and the action of ciJ is by — 1. The Fi— irreducible 

components of R^^ are Xk ©Pfc with actions {^k x ^k) © {^k x ^k)- Moreover, 
these are the isotypic components. The action of 6 on each component is by 
—1. As shown in 1 ^, the action of can be written in the form x + iy = 
e*®(x + iy). Therefore, vr G also acts by —1, so (l,'a7,7r) G Fi x Z 2 x 
acts trivially. 

There is a homomorphism Fi X Z2 X —)■ Fi X dehned by 


(p, 1,0) ^ (p,0) 

(p,OT,0) H- (p,0 + tt) 


and the action factors through this homomorphism, such that hnally there 
is a Fi X action, modulo K = {{l,zu,7i)). As shown in the isotropy 
subgroups are generated by the isotropy subgroups of the Fi x action 
together with K. 

Since we know what the action of Fi on C is, the corresponding action of 
Fi X is (^, 6) = uj'^^e^^z. The only way for S to have a two-dimensional 

hxed-point subspace is if S is the kernel of the Fi x action. This kernel 
consists of the elements (.^^, 9) such that = 1; that is, AmxTr/N + 

0 = 0. Hence S = {(.^^, —Arrixir/N)}. 

This group is isomorphic to Fi. Its spatial part S fl Fi consists of those 
elements for which my is divisible by N, so that y is a multiple of N/d 
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where d = gcd{m, N). That is, S fl Ti = Fd- Therefore 0121 ) represents 
a discrete rotating wave with spatial symmetry Fd- The discrete rotating 
wave is translated diagonally in the squared lattice i.e. the discrete rotating 
wave corresponding to the lattice nodes {a, a) with 0 < a < iV — a, is 
the same. This diagonal displacement is determined by the simultaneous 
and orthogonal action of a and p. As a consequence, there will exist N 
different discrete rotating waves per lattice direction a or /?; the number 
of total different discrete rotating waves per lattice is still N, because the 
discrete rotating waves in different rows (columns) of the lattice are just 
cyclic permutations of the N— waves of one row (column). 

Figure 3 shows the time evolution of the variables Xa ,/3 for a, = 1, 2, 3, 
for the case = 3 of one single torus. It can be seen that for hxed a, 
Xa,i 3 are three distinct discrete rotating waves, separated by a constant phase 
shift oi (j) = as predicted by the above analysis. Moreover, when switch¬ 
ing from a to a + 1, variable Xa+i,i 3 is again, phase shifted with ^ with 
respect to Xa,i 3 - In other words, discrete rotating waves are translated ’diago¬ 
nally’ in the lattice, so that we have Ha,p{t) = [xa^p{t — k(f)), — k(f))] = 

dd [a+k,y+k){mod Af) (^) ) with 0 ^ /c ^ N 1. 

5. Multifrequencies in a network of two conpled tori, each with Fi 

symmetry 

Multifrequencies patterns of oscillations have been initially observed in 
0 , 0. & and j^, where one or more oscillators oscillate at different 
frequencies. Golubitsky et ah and Stewart et al. developed the group- 
theoretical works which altogether represent the theoretical tools allowing a 
systematic analysis of these ^namics in systems with symmetry. Based on 
this theory, Palacios et al. [HI, [l^ as well as Longhini et al. [l^, devel¬ 
oped arguments that explained the experimental observation of these pat¬ 
terns in electrical circuits. We apply here Golubitsky’s theory and Palacios’s 
methodology to study the multifrequency patterns in a network obtained by 
diffusively coupling two tori. Goncretely, the Equation (|2]) becomes 
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(13) 


!/(«,/>),= Ko./J), - <5/(0,V’,1’=1,2 

Equation (fT^ represents the dynamical system of the two coupled tori; 
indices '0, v-, represent the two tori, e is the intertori coupling constant, while 
coupling functions i^{xa+i,j3,Xa,ii) and fi{xa,i3+i,Xa,i3) are dehned as in Equa¬ 
tion ([3]). 

Let us hrst observe that while each torus has interneuronal nearest- 
neighbor coupling in each direction, the intertori coupling is such that every 
neuron in one torus is coupled with a ’mean held’ of the neurons of 
the other torus. Therefore, while the symmetry group acting on each torus 
is X X Z 2 , the symmetry group of the full system is given by the 



2^ . which allows determining all oscillation patterns of the full system flT^ . 


The methodology consists in calculating the axial subgroups (up to conju- 



gacy) of the wreath product 


subgroups and maximal isotropy subgroups. It is a general approach able to 
theoretically predict the patterns of oscillation of a network formed by the 


nearest neighbor coupling of an arbitrary number of tori. This makes the 


object of another paper and is actually our work in progress that it will be 
published elsewhere shortly. In this Section, however, we do not pretend to 
carry out such a global analysis. The purpose of this section is to prove, 
by using a group theoretical argument, the possibility of a pattern of oscil¬ 
lations in which one torus produces discrete rotating waves with a constant 
phase shift, while the other torus shows in-phase oscillations at times 
the frequency of the rotating waves. We observed this multifrequency os¬ 
cillation pattern by computer simulations of system (IT^ and would like to 
theoretically prove the possibility of obtaining it. 

The idea is constructing a set of multifrequency patterns of interest and 
analyzing the action the symmetry group of the network has on the con- 
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structed set. It should be emphasized, the constructed set is only one of the 
possible sets of patterns of oscillations of the network. 

Let us hrst dehne 


Hi{i) 

H2{t) = [xaAt)^y<yAt)]T2 

to represent the state of the Xa,p,ya,j 3 variables in one tori 7^1 and 7 ^ 2 , re¬ 
spectively. Then the set P{t) = [iLi(t), Lf 2 (t)], represents an ensemble of the 
spatio-temporal patterns generated by the network. Computer simulations 
of system ([ 2 ]) corresponding to two coupled tori indicate that a possible set 
of patterns of oscillations is formed by discrete rotating waves with a con¬ 
stant phase shift in one torus and in-phase oscillations generated by the other 
torus. Figures |4](a) and|l](e), show three distinct traveling waves in two cou¬ 
pled tori of 3 X 3 neurons (Figure |4]( a)) and 11 x 11 neurons (Figure |4l(e)). 
The waves generated by torus are of the form 

HiTwit) = [xa,fiTw (^ - (« + A (mod N) 0), {t - {a + 0) (mod N) 0)] 

with a constant phase shift 0 = ^ -where r is the period of oscillations- 
among nearest-neighbor oscillators, the waves being cyclically permuted within 
the rows (columns) of the squared array. On the other hand, torus 702 gen¬ 
erates in-phase oscillations of the form 

\Xa,l3ip(t) T ya,l3 

of same period r, where all variables oscillate at N— times the frequency 
of the traveling waves and all [xa,/ 3 jp{t),ya,pjp(t)] are identical. A sample 
of the in-phase oscillations in two coupled tori of 3 x 3 neurons is shown in 
Figure 111(b) and for 11 x 11 neurons, in Figure l4](f). Therefore in this case, 
the set of patterns of oscillations P 2 Tif) = [Hirwif)) H 2 ip{t )], is formed by 
traveling waves {TW) and in-phase {IP) oscillations (where 2T stands for 
the case of two coupled tori). 

Now, it appears natural to search the spatial and temporal transforma¬ 
tions that leave unchanged the elements of the set of patterns of oscilla¬ 
tions. The set of these spatial and temporal transformations form the sym¬ 
metry group of the oscillation patterns [l^, [l^. As a hrst step, let’s as¬ 
sume that the pattern P 2 T = -^ 2 /p(^)] has the symmetry group 

(Fi X S^) X (Fi X S^), which describes (simultaneous) cyclic permutations 
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of the oscillators in each array, accompanied by shifts in time by 0. That is, 
(Fi X S^) ■ Hirpy^iit) = Hirpy^iit) so the traveling waves are unchanged. On 
the other hand, (Fi x §^) ■ H 2 ip{t) = H 2 ip{t + 0), so the in-phase oscilla¬ 
tors are shifted in time by (j). Thus if (Fi x S^) x (Fi x S^) is the symmetry 
group of P 2 T{t), then (Fi x x (Fi x §^) • P 2 t(^) = P 2 T{t), implies that 
H 2 ip{t) = H 2 jp{t + (p). There is a simple explanation for the fast oscillations 
generated by the second torus: the only way to obtain in-phase oscillations in 
torus #2 once we get traveling waves with constant phase shift in torus ^1, 
is when in-phase pattern oscillates at N times the frequency of the traveling 
wave pattern (l^ . 

When N is even, however, a condition for the Hopf bifurcation of the 
discrete rotating waves, taking into account is not only nearest-neighbor but 
also next-nearest neighbor couplings [l^,17|. As shown in jl^, there is a 


more subtle interpretation for bifurcation of discrete rotating waves when m 
and N are not coprime. 

As it was shown in Section H] and l3| , the spatial part S fl Fi consists 
of those elements for which my is divisible by N, so that y is a multiple of 
i = N/d where d = gcd{m, N). Then the ring of oscillators can be divided 
into L subsets of d oscillators, so that in each subset the oscillators behave 
identically (they oscillate in-phase), but the next subset oscillates with a 
phase-lag of T/l with respect to the hrst one. We therefore get a discrete 
rotating wave among sets of oscillators behaving identically. 


6. Concluding Remarks 

We analyzed the oscillatory dynamics of a network of electrically coupled 
FHN neurons. This network can be related to architectures of electrically 
coupled neurons observed in anatomical structures of the nervous system. We 
hrst described the building-block of the network, which consists of squared 
arrays shaped in the form of a torus. The analysis of the linear approximation 
of the system formed by one torus reveals the importance of interneuronal 
coupling strengths; they act as precursors of the Hopf bifurcation point. We 
give the analytical conditions for Hopf bifurcation and asymptotic stability 
of the coupled system, in which each torus is formed by an arbitrary number 
of neurons. We then performed a group theoretical analysis of one torus and 
found that only oscillatory pattern is represented by traveling waves with a 
constant phase shift. When analyzing two coupled tori, we identihed a set of 
patterns of oscillations whose elements are traveling waves and in-phase ’fast’ 
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oscillations. We showed that when the symmetry group of a single torus acts 
on the traveling waves, it leaves them unchanged; by the contrary, when it 
acts on the in-phase oscillations, they are shifted in time by (j). This proves 
that one of the possible patterns of oscillations of the network formed by 
the two coupled tori is represented by traveling waves produced by one torus 
and in-phase oscillations at iV— times the frequency of the traveling waves, 
shown by the other torus. This result is of possible interest for modeling the 
electrical activity of the nervous system. 
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Appendix A. Appendix 

Our aim is applying the criterium established by Guckenheimer and 
Holmes j^, to analyze the stability of the limit cycle of system ([T]) ob¬ 
tained by linearization about origin. From equation (E]), the condition for 
Hopf bifurcation is a = c< 6, a<l, and system o writes 


X = ax 
y = bx - 


- X 

ay. 


(A.l) 


It has been shown in [35 
be put in the form 


36| , that any system at Hopf bifurcation can 
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(A.2) 


X= -(py -I- f (x,y) 

y = (fx + g(x,y), 

by a suitable change of variables. In equation (1A.2p . is a real number, 
while / (x, y) and g {x, y) contain only higher-order nonlinear terms that 
vanish at origin. Moreover, one can discern between a subcritical and super¬ 
critical Hopf bifurcation by calculating the sign of the quantity: 


16s fxxx “1“ fxyy “h gxxy 9yyy 

1 (A.3) 

H [fxy ifXX A fyy') 9xy i9xx A 9yy) fxx9xx A fyy9yy\i 


where the subscripts indicate partial derivatives evaluated at (0,0). If 
s* < 0, the bifurcation is supercritical, while if s* > 0 the bifurcation is 
subcritical. 

In order to put system flA.lD in the form (IA.2H . we make the transforma¬ 
tion X = y - X, y = X. System flA.ip is transformed into 

V 6 — y/b — a? 


y = (px- f, 


a 


where (p = -\/b - a?, f {x, y) = g (x, y) = - f. Then, by eval- 

VO — CL 

uating expression flA.Op . we obtain 16s* = —6, which gives s* < 0 so Hopf 
bifurcation corresponding to steady state (0, 0), is supercritical. 
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two tori of two tori of 

3x3 neurons 3x3 neurons 




two tori of two tori of 


11 X 11 neurons 11 X 11 neurons 




two tori of 3 X 3 neurons each 


two tori of 11 X 11 neurons each 


Figure 4: Time series of variables (continuous line), 2:2,2 (dashed line) , 2:3,3 (dotted 
and dashed line) from the torus #1 (a) and torus #2 (b), from the integration of sys¬ 
tem m corresponding to two coupled tori, each of 3 x 3 neurons. Variables are phase 
shifted in torus #1, while they are synchronized in torus ^2. Coupling constants are, 
7 = 2 and S = 2, e = 0.5, while the rest of the parameters and initial conditions 
are the same as in Figure O Power Spectrum Density showing that neurons of torus 
#2 (d) oscillate at 1.29 = 3 x 0.431 Hz i.e. with 3-times the frequency of neurons 
from torus ^1 (c). Time series of variables 2 : 1,1 (continuous line), 2 : 5,6 (dashed line) 
, ;cii,ii (dotted and dashed line) from the torus #1 (e) and torus #2 (f), from the 
integration of system @ corresponding to two coupled tori, each of 11 x 11 neurons. 
Neurons are phase shifted in torus ^1, while they are synchronized in torus ^2. Cou¬ 
pling constants are, 7 = 8 and i5 = 8 , e = 0.5, while the rest of the parameters are 
the same as in Figure [2] Initial conditions of the neurons 2 :i„a, 2 /i„s, /? = 1,... ,4 of the 
torus #1 are xi,/j = (6.489,9.3843,6.9745,3.3656) = (0.8862,1.7536,8.3197,7.9935); 

initial conditions of the neurons 2 : 11 ,^ 3 ,j/n,^, (3 = 8 ,..., 11 of the torus #2 are 2 :i„a = 
(0.5475,9.7797,2.3655,3.8290) yij =’(4.8460,4.3027,7.5896,4.9924) while the initial’con¬ 
ditions for the rest of neurons in both tori have have been obtained by successive repetitions 
of the initial conditions used in Figure [2j Power Spectrum Density showing that neurons 
of torus 2 (h) oscillate at 57.1 = 11 x 0.519 Hz i.e. with 11-times the frequency of neurons 
from torus 1 (g). 
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